(Ga,Mn)As based superlattices and the search for antiferromagnetic interlayer 

coupling 
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Antiferromagnetic interlayer coupling in dilute magnetic semiconductor superlattices could result 
in the realisation of large magnetoresistance effects analogous to the giant magnetoresistance seen in 
metallic multilayer structures. In this paper we use a mean-field theory of carrier induced ferromag- 
netism to explore the multidimensional parameter space available in (Ga,Mn)As based superlattice 
systems. Based on these investigations we examine the feasibility of creating a superlattice that 
exhibits antiferromagnetic coupling and suggest potentially viable recipes. 
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I. INTRODUCTION 

The exciting new prospect of spin based electronics, 
known as spintronics, was initiated in 1988 with the 
discovery of giant magnetoresistance (GMR) in metallic 
multilayer structures These structures consist of in- 
terposed ferromagnetic (FM) and non-FM layers. When 
the magnetisation of adjacent FM layers is aligned in an- 
tiparallel directions, enhanced spin scattering of carriers 
causes an increased electrical resistance through the lay- 
ers, while when they are parallel the resistance is lower. 
Although typical GMR devices today consist of a tri- 
layer structure with a pinned magnetic layer and one in 
which the magnetisation is free to rotate, another method 
of implementation is with a superlattice structure where 
adjacent layers have an antiparallel magnetisation unless 
an external field is applied to align them. 

In multilayer structures containing ferromagnetic lay- 
ers, in addition to the ferromagnetic order within the 
layers, there can also exist magnetic exchange between 
the layers. The mechanism which causes the magnetic 
order between the layers is known as interlayer exchange 
coupling (IEC), and has been shown in metallic systems 
to be due to the spin polarisation of conduction carriers, 4 

Because the IEC energy considers the spin depen- 
dent changes in total energy, it thus determines which 
magnetic alignment of adjacent layers is energetically 
favourable. Although complicated helical arrangements 
can exist,— typically the interlayer exchange coupling will 
either be FM, where there is a parallel alignment of mag- 
netisation, or antiferromagnetic (AFM) where there is 
an antiparallel alignment. Therefore, in such a system, 
achieving AFM interlayer coupling is of high importance 
for technological applications. 

In addition to existing in metal systems, IEC 
is a generic property of magnetic multilayers, and 
AFM IEC has even been demonstrated in non-metallic 
FM semiconductor systems based on all-semiconductor 
EuS/PbS superlattices^ AFM IEC in dilute magnetic 
semiconductor (DMS) based superlattices was theoret- 
ically predicted in 1999 using a k • p kinetic- exchange 
model for carrier mediated ferromagnetism. 7 This ap- 



proach considers delocalised charge and adds extra mod- 
ulation induced by spin-polarised effects. A large 
magnetoresistance (MR) was predicted due to the large 
difference in miniband dispersion for the cases of fer- 
romagnetically and antiferromagnctically aligned lay- 
ers. Recently, IEC has been further explored using a 
tight-binding model.— This complementary microscopic 
approach, although not self-consistent, takes into ac- 
count atomic orbitals for all the constituent atoms, lead- 
ing to more accurate descriptions of the band struc- 
ture. Despite the different approaches used, both meth- 
ods provide qualitatively similar results for the IEC 
which shows oscillatory Ruderman-Kittel-Kasuya-Yosida 
(RKKY)-like behaviour. 

Although IEC has been shown to exist in DMS sys- 
tems based on (Ga,Mn)As/(Al,Ga)As trilayers^ there 
have been no reports of AFM interlayer coupling. Ex- 
perimental work into (Ga,Mn)As based multilayer and 
superlattice structures has only succeeded in demonstrat- 
ing FM IEC i 10 ^ 11 In order to test the prediction of a phe- 
nomenon analogous to GMR in metals in DMS materials, 
with a potentially much greater MR ratio, it is essential 
that AFM interlayer coupling is obtained. 

The aim of this study is to provide a comprehensive de- 
scription of the multidimensional parameter space avail- 
able in these DMS superlattice systems, in order to iden- 
tify optimal parameters for realising an antiferromagnet- 
ically coupled system. Because the interlayer coupling is 
mediated by carriers, a k ■ p approach is more practical 
for exploring a wide range of parameter values. The lim- 
itation of this approach is that a single parabolic band 
approximation is used, sacrificing full quantitative accu- 
racy for qualitative descriptions of a wide range of sys- 
tems. Subtleties of the band-structure and spin-orbit ef- 
fects are neglected. However, qualitative agreement with 
the data published in Ref. [8| at least partially justifies 
this approach. 

The organisation of this paper is as follows: first the 
details of the theoretical modelling of a DMS based su- 
perlattice system and the numerics of the self-consistent 
mean field calculations will be shown. Next, the re- 
sults, which will primarily consider (Ga,Mn)As based su- 
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perlattice systems with either GaAs or (Al,Ga)As non- 
magnetic spacer layers, will be presented. Finally, in the 
discussion, suggestions for recipes for superlattice sys- 
tems in which antiferromagnetic interlayer coupling may 
occur will be given. 
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which we shall rewrite as 



II. THEORETICAL MODELLING 

Our calculations are based on the Zener kinetic- 
exchange model 1 - description of magnetic interactions in 
Mn-doped III-V semiconductor structures. Microscop- 
ically, the kinetic- exchange between the local Mn mo- 
ments and itinerant hole spins originates from the p-d 
orbital hybridisation. 13 This model provides a good de- 
scription of ferromagnetism in bulk (Ga,Mn)As when the 
detailed structure of the valence band is taken into ac- 
count. 

An intuitive picture of the IEC in (III,Mn)V/III-V 
multilayer structures can be obtained by the perturba- 
tive mapping of the kinetic- exchange model onto an ef- 
fective interaction between local moments, following the 
RKKY approach, 7 The RKKY theory can be expected to 
provide useful predictions for structures close to a model 
pseudo-ID system consisting of alternating thin ferro- 
magnetic layers and non-magnetic spacer layers such that 
there is small coupling and low carrier polarisation. 14 
The RKKY range function falls off asymptotically with 
sm(2fc F d) ^ w j iere ^ j g carr j er wa ve vector and d is the 
distance between the magnetic layers. Thus, the RKKY 
theory shows that the coupling can have an oscillatory 
form. 

The Zener kinetic-exchange model for homogeneous 
(Ga,Mn)As was generalized in Ref. in order 
to account for the RKKY-like oscillatory effects in 
the inter-(Ga,Mn)As coupling in (Ga,Mn)As based 
ferromagnetic/non-magnetic superlattices on a more 
quantitative level. In this model the band structure is 
solved using the kinetic-exchange model and a parabolic 
band k ■ p effective mass approximation. In the Hamil- 
tonian the magnetic moments are accounted for through 
the p-d kinetic-exchange interaction between Mn spins 
and hole spins which is parametrized by a constant 
J p d and treated in the mean-held virtual crystal ap- 
proximation. The value of J pc i can be experimentally 
determined, and modern estimates of this value place 
it at 55 meVnm 3 £ To account for the inhomogene- 
ity, a standard formalisation of the local-spin density 
approximation (LSDA) using the Kohn-Sham equations 
for inhomogeneous systems is used in the band structure 
calculations, 1 ^ Hole mass is m* = 0.5m e and the spin of 
local Mn moments isS , = |atT = 0K. Thermodynam- 
ics are treated on a mean held level. 

In order to hnd the normalized wavefunction for 
a given energy, Bloch's theorem is used to solve 
the one-dimensional time-independent spin-dependent 
Schrodinger equation: 



*"(z) = /¥(*), (2) 

where k is the wavevector, n is the subband index, and 
a is the spin index and / = ^jp-(V a (z) — E). 
The Bloch function 

=u fc , n , ff (z)e< ifc *>, (3) 

gives the solutions of the Schrodinger equation for a 
periodic potential. 

For this system, the explicit form of the Hamiltonian 
for the spin-dependant potential V a {z) is given by 

V„{z) = V H + V xc<a + V b - a -\g*[i B B + h pd (z)}, (4) 

where Vjj is the Hartree (electrostatic) potential, given 
by the Poisson equation, V xc ,<x is the spin-dependent 
exchange-correlation potential given by the LSDA equa- 
tion, Vb is the band-offset, g* is the free-carrier g-factor 
and h p d is the mean-held kinetic-exchange interaction.— 
The IEC energy, E c , is defined as the difference in energy 
between the FM and AFM states per superlattice period. 

Let us suppose that our one-dimensional lattice has a 
period d n + m and consider now the solution only at N 
evenly distributed discrete points on the z-axis with a 
separation h. The wavefunction at each point is denoted 
as *&(z). By Taylor's theorem, the second-order approx- 
imations for \&(z + h) and ^(z — h) are 

V(z + h) = + M-'(z) + y 1<"(z); (5) 

9(z-h) = y(z)-hV'{z) + —y"(z). (6) 
Taking the sum of Eqs. (J3J) and © we obtain 

h 2 ^"{z) = #0 + h) + #0 -h)- 2*0). (7) 

Substituting the Schrodinger equation from Eq. ([2"|1 
into Eq. and rearranging gives the wavefunction at 
a given point as a linear combination of the wavefunc- 
tions at the two previous points: 

V(z + h) = (h 2 f + 2)y(z)-y(z-h). (8) 

This linear transformation can be represented as a 
transfer matrix, M„, such that 
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where ^ n is the wavefunction at the n z-point. By 
inspection we see that m n ,xx = h 2 f + 2, m n! i2 = — 1, 
™»i.2i = 1 and m nr 22 — 0. It is worth noting here that 
the determinant of each M„, det(M n ) = 1. The product 
of the N transfer matrices n , M n = M T represents 
the transformation from to ^n- By Bloch theorem's 
periodic boundary condition, Eq. ([3]) , this transformation 
can be written 




= e* 1 ( J; ) .(10) 

Therefore, 

= det(M T - e ika l 2 ) 

= (m?i - e lka ){m T 22 - e lka ) - m T l2 m T 2l 

= m\ x ml 2 - ml 2 ml x - e ika {m\ x + m T 22 ) + e 2lka 

= det(M T ) -e ifca Tr(M T ) + e 2ifca . (11) 

Since the determinant of M„ is 1, then the determi- 
nant of any product of M„, for any n, will also have a 
determinant of 1, hence dct(M T ) = 1. Substituting this 
into Eq. (|11[) gives 

1 - e ifca Tr(M T ) + e 2lka = (12) 

Tr(M T ) = e lfcQ +e~ jfca 

= 2cos(fca) (13) 

For a given energy, wavevector k can thus be found by 

k=- arccos(iTr(M T )), (14) 
a 2 

and the corresponding wavefunction can be found sim- 
ilarly. 



III. RESULTS 

In the RKKY model of interlayer exchange the oscilla- 
tions occur as a function of kp d, where kp is the Fermi 
wave vector and d is the separation between the two- 
dimensional magnetic planes^ 17 In our model we shall 



denote d n as the width of the non-magnetic layers, cor- 
responding to d from the RKKY model, and d m as the 
width of the magnetic layers. The length of a GaAs unit 
cell is labelled eto and has a value of 0.565 nm. We shall 
also define the average Fermi wave vector kp as 

k F = (3n 2 N 3D )i, (15) 

corresponding to the Fermi vector kp in the ideal 
RKKY model with a parabolic band. The average 3D 
carrier concentration N 3 p> is defined as 

N 3D = — !— / N 3D (z) dz = j^. (16) 

"n+m J unit coll O-n+rn 

The superlattice structures being considered in this pa- 
per consist of thin (Ga,Mn)As layers interposed with non- 
magnetic spacer layers. The primary structural parame- 
ters that can be changed are the widths of the layers and 
their composition. Fig. [Ija) shows the calculated self- 
consistent charge distribution and potentials for a sim- 
ple case with a low moment concentration (2%), and the 
spacer layers are thicker than the magnetic layers. There 
is a uniform impurity concentration of acceptors through- 
out the structure, either magnetic, like Mn in the mag- 
netic layers, or from non-magnetic dopants in the spacer 
layers. The polarisation of carriers is low in this case, and 
they have an almost totally uniform distribution. These 
figures show cases where the number of monolayers of 
non-magnetic layer d n j \a§ = 5 and average 3D carrier 
concentration N$p> — 10 20 cm" 3 . The Fermi energy is 
at V a — eV. The interlayer coupling is in an assumed 
AFM state, although this may not be the energetically 
favoured state for such a system. 

The effect of increasing the number of moments in the 
magnetic layers is shown in Fig. [l^b) , where the Mn dop- 
ing is 8%. As expected, the increased moment concen- 
tration increases the size of spin splitting in the magnetic 
layers, which increases the carrier polarisation. Addition- 
ally, the potential due to the magnetic ordering causes a 
redistribution of carriers to occur, increasing the concen- 
tration in the magnetic layer. Another way to cause car- 
rier redistribution is to remove the doping from the non- 
magnetic layer, which is shown in Fig.[TJc), where the Mn 
concentration in the magnetic layers is again 2%. With- 
out a neutralising background charge, Coulomb repulsion 
opposes carrier redistribution into the non-magnetic lay- 
ers. The resulting carrier distribution is similar to that 
of Fig. [IJb). Although there is a greater concentration 
of carriers in the magnetic layers, the polarisation is not 
significantly increased over that seen in Fig. [IJa) where 
there is a uniform charge concentration. 

To cause stronger confinement of carrier to the mag- 
netic layers, the non-magnetic layers can be made from 
(Al,Ga)As. The effects of this are shown in Fig. Ufa), 
where the carrier concentration in the centre of the mag- 
netic layers is about double that of the centre of the non- 
magnetic layers. Although there is a slight increase in 
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FIG. 1: (colour online) The self-consistent charge distribution, N 3 £>, and potentials, V a , for a double unit cell of 
three different (Ga,Mn)As/GaAs based superlattice structures in an AFM state. d m /^a — 2, d n /^a — 5 and 

iV 3 £> = 10 20 cm -3 in each case. V a = eV corresponds to the Fermi level, (a) 2% Mn doping and a uniform 
impurity concentration, (b) 8% Mn doping and a uniform impurity concentration and (c) 2% Mn doping but no 

impurities in the non-magnetic layer. 



the polarisation of carriers over the previous two cases 
with 2% Mn moment, this is increased significantly when 
more moments (8%) are put in the magnetic layers, as 
shown in Fig. El^b). Because of the strong confinement 
of carriers to the highly magnetic layers there is a very 
high polarisation of carriers; this effect is enhanced over 
that of the 8% case with a doped GaAs spacer which 
was shown in Fig. [T^b) . Finally, Fig. &[c) shows a case 
where the non-magnetic layers are thinner than the mag- 
netic layers, and also have a high Mn doping. Following 
the established trend, carriers are strongly polarised and 
tightly confined to the magnetic layers. In fact, with an 
almost total polarisation and a very high depletion of the 
non-magnetic layers, this structure represents an almost 
opposite case to that of Fig. [lja) . 



Bearing in mind these examples of how changing the 
structural properties can alter the electronic configura- 
tion of the superlattices, the effects of these changes on 
the IEC will now be explored. This will be done in 
two parts. Firstly, GaAs based spacers, similar to those 
shown in Fig. [T] will be considered. The more extreme 
cases presented in Fig. [2] will be considered in the second 
half. 



A. GaAs spacer 

First to be considered is a superlattice structure close 
to the RKKY limit of infinitely thin magnetic layers sur- 
rounded by free unpolarised carriers. For this we shall 
use thin magnetic layers and a low magnetic moment con- 
centration, as per Fig.QJa). In Fig.[3^a) the IEC energy, 
E c , is plotted against the 3D carrier concentration, N3D, 
and number of monolayers of GaAs in the non-magnetic 
spacer, 2d n /^ao- The magnetic (Ga,Mn)As layer is 2 
monolayers thick and contains 2% Mn local moment dop- 
ing. There is a uniform acceptor density throughout the 
structure which gives an average hole concentration of 
4.43 x 10 20 cm~ 3 . In this case there are oscillations as a 
function of both parameters, analogous to the kpd oscil- 
lations in the ideal quasi one-dimensional RKKY model. 
For the calculated IEC energy, E c , positive values cor- 
respond to FM interlayer coupling being energetically 
favourable, and negative values correspond to AFM in- 
terlayer coupling being the favoured configuration. 

The RKKY like behaviour observed in Fig. [2^a) is 
consistent with the results obtained in the tight-binding 
approach^ when the exchange coupling energy, E c , is 
plotted against the two-dimensional carrier concentra- 
tion, N2D1 for fixed layer thicknesses. However, it is 
worth noting that when the exchange coupling is plot- 
ted as a function of the non-magnetic spacer thickness, 
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FIG. 2: (colour online) The self-consistent charge distribution, N 3 d, and potentials, V a , for a double unit cell of 
three different (Ga,Mn)As/(Al,Ga)As based superlattice structures in an AFM state. d n /^ao — 5, N^d — 10 20 cm - ' 
and the Al concentration is 30% in each case. V a = eV corresponds to the Fermi level, (a) d m /\a$ = 2 and 2% 
Mn doping (b) d m /^ao = 2 and 8% Mn doping (c) d m /^ao = 8 and 8% Mn doping. 



d n , for a fixed N2D there are no apparent RKKY oscilla- 
tions. Because N^d, and therefore kp, is a function of d n 
these two parameters are not independent when A^d is 
fixed. This results in the oscillatory behaviour appearing 
to be suppressed. 

There are, however, real physical reasons for deviation 
from RKKY behaviour. The data from Fig. OUa) are 
rcplotted in Fig.[5]as a function of 2kpd n +i- Also plotted 
is the function 



where a is a scaling factor. This function is the asymp- 
totic limit of the pseudo one-dimensional RKKY range 
function^ The strength of the interaction is expected 
to scale with the density of states, and in the ID case 
a ~ k F i^ The different series of points on the graph 
correspond to the series of different non-magnetic spacer 
thicknesses shown in Fig. [3]Ja). For a given 2kpd n +i, 
the points with the largest magnitude are those with 
the greatest /c_f; this behaviour is consistent with the 
expected scaling of a with kp. Note the fact that, in 
order to have improved alignment of the curves, the os- 
cillations were plotted with the effective value of the non- 
magnetic spacer, d ni being increased in size by one mono- 
layer, which is denoted d n +\. This is necessary due to 
the fact that the magnetic layer in this structure is not 



infinitely thin, as per the ideal RKKY case, but has a 
defined width. 

Exploring this deviation from RKKY behaviour fur- 
ther, Fig. |3Kb) shows the IEC for a superlattice system 
with a thicker magnetic layer, d m /\aQ = 8. All other pa- 
rameters are as with (a). Examining the AFM peak, the 
reduction in average carrier concentration of the mini- 
mum as the spacer thickness is increased occurs more 
rapidly, evidenced by the larger derivative of the aver- 
age 3D carrier density, N^p>, with respect to the non- 
magnetic layer thinkness, d n , of the minimum E c at low 
d n . By way of contrast, at large d n this is lower, that is 
the curve has become much more straight. This is con- 
sistent with the effects of large magnetic layers increas- 
ing the centre-to-centre distance of the magnetic layers, 
causing the effect of an apparently larger non-magnetic 
layer. However, in addition to this, increasing the mag- 
netic layer thickness has introduced additional points of 
inflection, for reasons that are not immediately obvious. 

It is also possible to deviate from RKKY-type be- 
haviour through redistribution of charge. There are two 
primary methods by which this is achieved, as shown in 
Figs. [lib) and (c). The first is that charge is confined 
to the magnetic layers by the magnetic exchange poten- 
tial. Fig. HJa) shows the IEC where the Mn doping has 
been increased to 8%. However, when the magnetic layer 
is thin, significant charge redistribution is opposed by 
the Coulomb potential and the RKKY character is not 
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(a) d m /iao = 2, 8% Mn concentration 
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(b) d m /^ao = 8, 2% Mn concentration 

FIG. 3: (colour online) The IEC energy, E c , of two 
(Ga,Mn)As/GaAs based superlattices as a function of 

the average 3D carrier concentration, N^d, and the 
number of monolayers of non-magnetic layer, 2d n j\aQ. 
Positive (red coloured) values of E c indicate FM 
intcrlayer coupling is energetically favourable and 
negative (blue coloured) values indicate AFM is 
favourable. Both superlattices have a 2% Mn doping in 
the magnetic layer and there is a uniform impurity 
concentration of acceptors through the structure. 



significantly affected. As the figure shows, the main ef- 
fect is that the size of the IEC is increased. Despite the 
larger spin splitting and the larger polarisation of car- 
riers caused by the greater moment concentration, the 
coupling retains an RKKY character. 

When the magnetic layer is made wider the increased 
quantity of magnetic moments now causes additional 
changes in the oscillatory behaviour, beyond that of sim- 
ply increasing d m . Fig. 0|b) plots the IEC for a system 
which now has magnetic layers of 8 monolayers with a 
Mn doping of 8%. Because of the increased depletion of 
carriers from the non-magnetic layers, the N30 values at 



(b) d m /^ao = 8, 8% Mn concentration 

FIG. 4: (colour online) The IEC energy, E c , of two 
(Ga,Mn)As/GaAs based superlattices as a function of 

the average 3D carrier concentration, N3D, and the 
number of monolayers of non-magnetic layer, 2d n j\a§. 
Both superlattices have an 8% Mn doping in the 
magnetic layer and there is a uniform impurity 
concentration of acceptors through the structure. 



which AFM coupling is expected to occur are now greater 
for a given non-magnetic layer thickness. Additionally, 
the damping of the magnitude of the IEC energy with 
increasing d n has now significantly changed. While the 
first FM and AFM maxima are rapidly diminished with 
increasing non-magnetic spacer, the second FM peak is 
not greatly affected. The second AFM peak even in- 
creases in magnitude with larger d n , and for large spacer 
it can even be greater than the first. 

Note that when the unit cell becomes large and there 
is a high carrier concentration, the weak coupling and 
flat minibands make self-consistent convergence difficult; 
these regions are visible as rough areas on the figures. 
Such samples would anyway be extremely sensitive to 
inhomogeneities and fluctations. No data is shown where 
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2kpd n -\-i 

FIG. 5: (colour online) IEC energy, E c , as a function of 
2kpd n+ i for a superlattice with magnetic layers with a 

Mn doping of 2% and 2 monolayer thickness, and a 
uniform impurity concentration. The (red) curve is an 
estimate of the ideal RKKY range function. 



the calculations have diverged. 

The second method of charge redistribution is via a 
Coulomb potential. Fig. [D(a) shows the IEC profile for a 
system with a magnetic spacer of two monolayers and a 
Mn concentration of 2%. However, now there is no neu- 
tralising background charge in the non-magnetic layer, 
so self-consistent redistribution results in the formation 
of an effective barrier. Fig. QJc) shows the potentials 
and charge distribution for a unit cell of this structure 
in an AFM configuration, again with d n /^ao = 5 and 
N^E) = 10 20 cm . As previously noted, the Coulomb 
barrier formed is comparable in size to the spin splitting 
caused by the 2% Mn doping. This results in a similar 
charge redistribution as in the 8% doped case, although 
without such strong carrier polarization. As with that 
case, there is not a significant deviation from RKKY- 
type behaviour. 

Increasing the magnetic spacer thickness now causes 
more significant changes than seen with the doped spac- 
ers. Fig. [U(b) shows the IEC profile for a superlattice 
with d m /^aQ = 8 with a 2% Mn doping and no impuri- 
ties in the non-magnetic spacer. In addition to the extra 
inflection points there is now an additional AFM region. 
The magnitude of the local minimum in this region does 
not decrease much with non-magnetic spacer width, and 
occurs with an almost linear dN$D /dd n ■ This is now 
very unlike RKKY behaviour. 

To investigate this further we shall now consider su- 
perlattice with (Al,Ga)As non-magnetic spacers, so that 
greater charge redistribution will occur than that caused 
by the magnetic ordering potential of a high magnetic 
moment concentration, or the Coulomb potential arising 
from an undoped spacer. 
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(a) dm/h a = 2, 2% Mn concentration 
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(b) d m /^ao = 8, 2% Mn concentration 

FIG. 6: (colour online) The IEC energy, E c of two 
(Ga,Mn)As/GaAs based superlattices as a function of 

the average 3D carrier concentration, N^d, and the 
number of monolayers of non-magnetic layer, 2d n /^ao. 
Both superlattices have a 2% Mn doping in the 
magnetic layer but there is no charge doping in the 
non-magnetic spacer layer. 



B. (Al,Ga)As spacer 

In the previous section it was demonstrated that inter- 
layer coupling in superlattice structures would have an 
oscillatory behaviour as a function of parameters N^jj 
and d n , analogous to that of RKKY, when the magnetic 
layers were thin and surrounded by charge. As the struc- 
ture of the superlattice is changed the IEC would start to 
deviate from the ideal RKKY behaviour. This is partic- 
ularly apparent with increased magnetic layer thickness. 
Changing the 3D charge distribution has a more limited 
effect; neither large magnetic moment concentration nor 
a self-consistent Coulomb barrier would cause significant 
confinement of carriers. In order to investigate these 
effects further, a band offset will be introduced to fur- 
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(a) rfm/ido = 2, 2% Mn concentration 



(a) d m /iao = 2, 8% Mn concentration 
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(b) dml\a,Q = 8, 2% Mn concentration 

FIG. 7: (colour online) The IEC, E c , of two 
(Ga,Mn)As/(Al,Ga)As based superlattices as a function 
of the average 3D carrier concentration, N 3 £>, and the 
number of monolayers of non-magnetic layer, 2d n /^ao. 
Both superlattices have a 2% Mn doping in the 
magnetic layer and the non-magnetic layers are 
(Al .3,Gao.7)As and have no charge doping. 



(b) dm/^ao = 8, 8% Mn concentration 

FIG. 8: (colour online) The IEC, E c , of two 
(Ga,Mn)As/(Al,Ga)As based superlattices as a function 
of the average 3D carrier concentration, N313, and the 
number of monolayers of non-magnetic layer, 2d n /^ao. 
Both superlattices have an 8% Mn doping in the 
magnetic layer and the non-magnetic layers are 
(Al .3,Gao.7)As and have no charge doping. 



ther confine carriers to the magnetic layers. This will be 
achieved by using (Alo.3,Gao.7)As as the non-magnetic 
layer material, which has a valence band offset of about 
150 meV from GaAsJ^ 

Fig. EKa) shows the IEC profile for a structure with 
a (Gao.98,Mno.o2)As magnetic layer of 2 monolayers and 
an (Alo.3,Gao.7)As non-magnetic layer. There is no dop- 
ing in the non-magnetic layers. The peak FM and AFM 
coupling strengths are now stronger than in the case with 
doped GaAs spacers seen in the otherwise identical struc- 
ture in Fig.[3](a). The charge distribution of this structure 
is shown in Fig.[2fa), where the barrier confines carriers 
to the magnetic layers, as expected. However, the 2d n kp 
oscillations are damped more rapidly than with the GaAs 
spacer, resulting in the second FM and AFM peaks being 



very weak. This additional damping occurs particularly 
rapidly with increasing carrier density, N^o- As a result, 
the first anti-ferromagnetic peak barely reduces in mag- 
nitude as the non-magnetic layer thickness is increased. 
Juxtaposing this with the GaAs barrier case, where the 
largest AFM coupling that can occur when d n /^a = 10 
is less than a quarter of the size of that when d n /^ao = 2, 
we identify this as a significant departure from the pre- 
viously seen RKKY-like oscillatory behaviour. 

Increasing the magnetic moment concentration leads 
to a more interesting alteration than with the GaAs 
spacer, where the effect was principally to scale up the 
magnitude of the IEC energy. Fig. [5Ja) shows the IEC 
for a (Ga,Mn)As/(Al,Ga)As superlattice with an 8% Mn 
doping in the 2 monolayer magnetic layer, as previously 
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considered in Fig.[2jb). Now the first AFM peak appears 
to have two stages. The first is at low spacer thicknesses, 
where the average hole density at which the maximum 
occurs decreases with increasing spacer thickness. For 
large spacer thicknesses the curve has straightened out, 
and there is almost no dependence on d n for the sign of 
the coupling. This characteristic is similar to that ex- 
hibited in Figs.[3Jb) and[Hb), where the magnetic layer 
is 8 monolayers thick. This was attributed to loss of 
independence of the d n and N$d parameters, as the sys- 
tem became less RKKY-like. Knowing that the band off- 
set and large magnetic ordering cause significant carrier 
redistribution, particularly, this means that the carrier 
concentration in the spacer will decrease as a function of 
spacer thickness. This can account for the weak depen- 
dence of E c on kpd n . Also, note that the size of the first 
AFM peak decreases more rapidly at high spacer thick- 
nesses where the average carrier concentration, Nsd, at 
which it occurs is not decreasing. This is consistent with 
the previous observation of enhanced damping with in- 
creasing carrier concentration. 

With high magnetic layer thicknesses the RKKY-type 
oscillations have almost completely disappeared. The 
beating patterns which were emerging in the d m /\aa = 8 
GaAs spacer cases have now come to dominate the IEC. 
Figs.[7jb) and[8]Jb) shows this for d m j\aQ = 8, with re- 
spectively 2% and 8% Mn doping. In these cases the 
oscillations occur almost exclusively as a function of hole 
density, being almost independent of the spacer thick- 
ness. Note, however, as was seen in Fig. [2jc), the non- 
magnetic layer is highly depleted when the magnetic layer 
is 8 monolayers thick with an 8% Mn doping. This makes 
computing IEC for larger spacers unfeasible. 



IV. DISCUSSION AND RECIPES 

Having explored the parameter spaces we will now con- 
sider possible structures of a (Ga,Mn)As based superlat- 
tice that would exhibit AFM interlayer coupling. Each 
parameter will be considered for feasibility, and, based 
on the above calculations, suggestions for values can be 
made. 

The first to be considered is the Mn concentration in 
the (Ga,Mn)As layers. From the viewpoint of simply cre- 
ating a viable ferromagnet this is an essential parameter; 
not only does each substitutional Mn provide a magnetic 
moment, it also acts as an acceptor and thus this fac- 
tor controls the hole concentration. Calculations 2 ^ esti- 
mate that the minimum hole density for ferromagnetism 
is ~ I0 20 cm~ 3 . Assuming that each Mn provides one 
hole, this carrier concentration would correspond to a 
moment concentration of ~ 0.5%. Experimentally, typi- 
cal Mn concentrations are in the range of 2-8% (4.4 x 10 20 
to 1.8 x 10 21 cm -3 ), which in good quality material could 
result in higher carrier concentrations than the 10 to 
10 21 cm -3 range considered in these calculations. Whilst 
the higher magnetic moment concentration can increase 



the size of the E c peak, and thus a high moment concen- 
tration is favourable, the high carrier concentrations that 
would be associated with this would cause the strength 
of the IEC to become extremely weak. This constraint 
therefore imposes a practical range for Mn concentrations 
as being between 2 and 4% (4.4 x 10 20 to 8.8 x 10 20 cm -3 
respectively) . 

For the non-magnetic spacer thickness the general 
trend is that the strength of the IEC becomes weaker as 
the non-magnetic layer becomes thicker. Although this 
effect is somewhat diminished for the cases where there 
is strong carrier confinement to the magnetic layers, it 
is a serious consideration and, ideally, to see strong IEC 
effects, the non-magnetic layer should be as thin as pos- 
sible. Furthermore, particularly in cases where the 2kpd 
behaviour is dominant, as carrier concentration increases 
the spacer thickness at which the AFM IEC is strongest 
decreases inversely. As discussed above, low carrier con- 
centrations are not possible, so therefore it would seem 
beneficial to make the spacer layers as thin as practicable. 
Bearing in mind that the average distance between two 
Mn atoms when the concentration is 3% is of the order 
of a couple of GaAs unit cells, in order to make the non- 
magnetic spacer a discernible barrier then 4 monolayers 
would seem to be a realistic lower bound. 

The effect of the magnetic layer thickness on the IEC 
profile is more subtle, and seems mainly to distort the 
RKKY behaviour but otherwise, in the limits consid- 
ered within this study, does not have any negative ef- 
fects on the interlayer coupling. However, again for in- 
terlayer coupling to exist it is necessary that each mag- 
netic layer is itself ferromagnetic. Usually (Ga,Mn)As is 
grown in bulk layers of many nanometres; the thinnest 
(Ga,Mn)As epilayers for which published literature exists 
are 5 nm thick; 22 and the (Ga,Mn)As based heterostruc- 
tures with magnetic layers as thin as 8 monolayers have 
been reported to be ferromagnetic. 23 It would therefore 
seem prudent, in order to ensure that the magnetic lay- 
ers are effective ferromagnets, to prefer to make them 
thicker. For the 5 nm film some amount of surface deple- 
tion should be expected, so a 8 monolayer thick magnetic 
layer, equivalent to 2.26 nm, is comparable. Of course, 
if thinner films are shown to be viable then there is no 
reason not to consider them also. 

Based on these constraints, Fig. [5] shows the IEC pro- 
file for two candidate superlattices as a function of car- 
rier concentration, N$d- Both superlattices are identi- 
cal in structure except for the composition of the non- 
magnetic layer. The magnetic layer thickness is 8 mono- 
layers and has a magnetic impurity concentration of 
5 x 10 20 cm -3 (2.26%) and the non-magnetic layers are 
4 monolayers thick. As expected from the calculations, 
when the (Al,Ga)As barriers strongly confine carriers 
to the magnetic layers the IEC energy can have poten- 
tially greater magnitudes, although the oscillations have 
a much higher frequency. In these samples the carrier 
concentration would be somewhere below the Mn con- 
centration of 5 x 10 20 cm -3 (2.26%), however the exact 
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FIG. 9: (colour online) A comparison of the IEC energy, 

J? c , as a function of the average 3D carrier 
concentration, N3D, for two specific superlattices with 
either a GaAs or an (Alo.3,Gao.7)As non-magnetic layer. 
The magnetic layers are 8 monolayers thick and have a 
Mn concentration of 5 x 10 20 cm -3 (2.26%) and the 
non-magnetic layer is 4 monolayers thick. 



amount would depend on subtleties of the growth condi- 
tions. Although this suggests that for AFM IEC to occur 
the desired carrier concentration should be several times 
lower, it must be accepted that the calculations are of a 
more qualitative nature. Additionally, by tailoring the 
band offset of the non-magnetic layer by altering the Al 
content, the location of the peak can be adjusted some- 
what. This at least shows that these designs offer the 
possibility for AFM interlayer coupling. 

Even if the IEC energy were to favour an AFM ar- 
rangement, if the AFM coupling is weaker than the 
anisotropy fields it is possible that, after the applica- 
tion of a field, the superlattice could become locked into 
a FM spin configuration. This spin-locking behaviour 
has been observed in EuS/PbS superlattices^ and Fe/Nb 
multilayers 2 ^ studied via neutron scattering. 

Comparing, then, the calculated IEC to the mag- 
netocrystalline anisotropic energy of (Ga,Mn)As, we 
take a typical "worst case" value of the in-plane cubic 
anisotropy constant to be of the order of 2000 J m -3 at 
4.2 K.— Using a value of the interlayer coupling energy 
E c = 10 uJm~ 2 from Fig. [5] and using the bilayer period 
of 3.4 nm we find the energy density of the IEC energy 
is 3000 J m~ 3 . Although this is assuming an ideal value 



of E c , this compares favourably with the anisotropy en- 
ergy. Furthermore, larger values for the IEC have been 
found in the tight-binding approach.— Therefore, such a 
superlattice structure might reasonably be expected to 
be a candidate to exhibit AFM interlayer coupling. 
V. CONCLUSION 

The composition and structure of (Ga,Mn)As based 
superlattices can have profound effects on the expected 
IEC. By examining possible compositions within the 
broad parameter space that these structures offer it is 
possible to identify different recipes for devices that 
might offer the possibility of demonstrating AFM in- 
terlayer coupling. Such a study was presented in this 
paper based on the parabolic band k • p kinetic ex- 
change model. This model ignores spin-orbit and band 
warping effects, although comparisons to the previously 
studied microscopic tight-binding model suggest that 
our results provide a reasonable qualitative or semi- 
quantitative description of the system. The calculations 
predict that short period superlattices with magnetic and 
non-magnetic layers with widths less than 10 monolay- 
ers seem be promising candidates. There is existing ex- 
perimental work in (Ga,Mn)As/GaAs based superlattices 
with similar dimensions, but this has only exhibited FM 
IE G 23 i 26 The ideal dimensions suggested by the calcula- 
tions are therefore feasible and have shown themselves to 
be viable ferromagnetic DMS material. 

Particularly interestingly, in a (Ga,Mn)As/(In,Ga)As 
based superlattice composed of 8 and 4 monolayer width 
layers respectively two phase transitions were observed^ 
Although there was no evidence that this was due to any 
AFM effects, this indicates that there is some additional 
physics at play in these systems, making them of interest 
for further study. The calculations presented in this pa- 
per have shown the importance of the composition of the 
non-magnetic spacer on the character of the IEC. Us- 
ing different Al concentrations to tailor the band-offset 
between layers in a superlattice is a standard tool for de- 
signing normal non-magnetic superlattice systems. Util- 
ising this technique in the magnetic superlattices could 
potentially provide a way to tune to IEC profile to one 
where AFM coupling is preferential. 
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